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Ion bombardment is known to enhance surface diffusion and affect the surface morphology. To 
quantify this phenomenon we calculate the ion-induced diffusion constant and its dependence on 
the ion energy, flux and angle of incidence. We find that ion bombardment can both enhance 
and suppress diffusion and that the sign of the diffusion constant depends on the experimental 
parameters. The effect of ion-induced diffusion on ripple formation and roughening of ion-sputtered 
surfaces is discussed and summarized in a morphological phase diagram. 

PACS numbers: 79.20.Rf, 64.60. Ht, 68.35.Rh 



Sputtering, the removal of atoms from the surface of 
solids through the impact of energetic particles (ions), is 
an important thin film processing technique Conse- 
quently, much attention has been focused on the mea- 
surement and calculation of the sputtering yield and of 
the velocity and angular distribution of the sputtered 
particles. However, for many applications an equally 
important phenomenon, ion-induced surface diffusion, 
has eluded sufficient understanding so far. In the ab- 
sence of ion bombardment surface diffusion is thermally 
activated and characterized by the diffusion constant 
Dt = Doexp[—Ed/kBT] such that the evolution of the 
surface height h(x,y,t) is described by the continuum 
equation dh/dt ~ — I>rV 4 /i Here Ed is the activation 
energy for surface diffusion of the adatoms and T is the 
substrate temperature. However, numerous experiments 
investigated the effect of ion bombardment on island for- 
mation, surface migration, surface smoothing and ripple 
formation have provided evidence that ion bombardment 
is accompanied by an increase in surface diffusion |3] |). 
In particular, it has been demonstrated that ion-induced 
surface diffusion can decrease the epitaxial temperature 
H , enhance nucleation during growth Q] , or modify the 
surface morphology. For example, MacLaren et al. |^] 
bombarded GaAs with 17 KeV Cs + in the temperature 
range from —50 to 200 °C, observing the development of a 
ripple structure on the surface with a wavelength propor- 
tional to the square root of the diffusion constant. When 
decreasing the temperature, the ripple spacing (wave- 
length) did not decrease exponentially to zero with the 
inverse temperature, but at approximately 60 °C it sta- 
bilized at a constant value, providing direct evidence for 
a temperature independent ion-induced surface diffusion 
constant. 

Although the effect of the ions on surface diffusion is 
well documented experimentally, there is no theory that 
would quantify it. In this paper we calculate analytically 
the ion- induced diffusion constant, D 1 , and its depen- 
dence on the ion energy, flux, angle of incidence, and 
penetration depth. We find that there exists a param- 
eter range when ion bombardment generates a negative 
surface diffusion constant, leading to morphological in- 
stabilities along the surface, affecting the surface rough- 



ness and the ripple structure. The effect of ion-induced 
diffusion on the morphology of ion-sputtered surfaces is 
summarized in a morphological phase diagram, allowing 
for direct experimental verification of our predictions. 

Ion-sputtering is determined by atomic processes tak- 
ing place within a finite penetration depth inside the 
bombarded material. The ions penetrate the surface and 
transfer their kinetic energy to the atoms of the substrate 
by inducing cascades of collisions among the substrate 
atoms, or through other processes, such as electronic ex- 
citations. A convenient picture of the ion bombardment 
process is shown in Fig. [j]. The ions penetrate a distance 
a inside the solid before they completely spread out their 
kinetic energy with some assumed spatial distribution. 
An ion releasing its energy at point P in the solid con- 
tributes energy to the surface point O, that may induce 
the atoms in O to break their bonds and either leave the 
surface or diffuse along it. 




FIG. 1. Following a straight trajectory (solid line) the ion 
penetrates an average distance a inside the solid (dotted line) 
after which it completely spreads out its kinetic energy. The 
energy decreases with the distance from P, the dotted curves 
indicating schematically the equal energy contours. The en- 
ergy released at point P contributes to erosion at O. The 
inset shows the laboratory coordinate frame: the ion beam 
forms an angle 6 with the normal to the average surface ori- 
entation, z, and the in-plane direction x is chosen along the 
projection of the ion beam. 
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Following ||,|l0|], we consider that the average energy 
deposited at point O due to the ion arriving at P follows 
the Gaussian distribution 
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In (|l|) z' is the distance measured along the ion trajec- 
tory, x' , y' are measured in the plane perpendicular to it; 
e denotes the kinetic energy of the ion and a and /x are the 
widths of the distribution in directions parallel and per- 
pendicular to the incoming beam, respectively. However, 
the sample is subject to an uniform flux J of bombarding 
ions. A large number of ions penetrate the solid at dif- 
ferent points simultaneously and the velocity of erosion 
at O depends on the total power So contributed by all 
the ions deposited within the range 1Z of the distribution 
(HI), such that 
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where $(r) corrects for the local slope dependence of the 
uniform flux J and p is a proportionality constant be- 
tween power deposition and the erosion rate ||. The 
calculation of v involves the following assumptions (for 
more details see Refs. [10 111 ): (a) In the laboratory co- 



ordinate frame (x, y, z) the surface can be described by a 
single valued height function h(x, y, t), measured from an 
initial flat configuration which lies in the (x,y) plane (see 
Fig. Q); (b) The angle between the ion beam direction 
and the local normal to the surface is a function of the 
angle of incidence 9 and the values of the local slopes d x h 
and d y h, and can be expanded in powers of the latter. 
Under these conditions, we can expand (Q), obtaining the 
equation of motion [fL2| 
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From (g) we can calculate the expressions for the co- 
efficients appearing in (|J) in terms of the physical pa- 
rameters which characterize the sputtering process. The 
coefficients v x and v y were first calculated by Bradley 
and Harpe r |10| , while the nonlinear expansion was per- 
formed in |llf|, providing \ x and X y . A fourth order ex- 
pansion is used to obtain the ion-induced diffusion con- 
stants D\. and Dy. To simplify the discussion we re- 
strict ourselves to the symmetric case a = p,. Using 
F = (eJp/V2~7r) exp(-a 2 /2 + a 2 s 2 ), s = sin#, c = cosO 
and a a = a/ a, we find for the surface diffusion constants 
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Consistent with symmetry considerations for 9 — we 
obtain D\. = Dy. However, for 9 we find that 
D[ D' y , i.e. the ion-induced surface diffusion is 
anisotropic. Moreover, its sign also depends on the ex- 
perimental parameters. The consequences of (^) and (|5|) 
can be summarized as follows: (a) Independent of the 
angle of incidence Dy is positive while the sign of the D^. 
depends on both 9 and a a as shown in Fig. @. 
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FIG. 2. Ion-induced diffusion constant, Di and Dy (inset) 
as a function of the angle of incidence 6. In both figures the 
curves correspond to a CT = 1.0 (dashed line), a a = 1.5 (dotted 
line) and a a = 2.0 (continuous line). 

Thus, while for 9 = the ion bombardment enhances 
the surface diffusion (D^, > 0), for large 9 it can suppress 
diffusion; (b) The fact that D\. can be negative indicates 
that any simple theory connecting the magnitude of the 
ion-induced diffusion to the energy transferred by the 
ions to the surface is incomplete, since it can predict only 
a positive D 1 . In fact, D 1 is the result of a complex in- 
terplay between the local surface topology and the energy 
transferred to the surface; (c) The diffusion constants are 
proportional to the flux J, in agreement with the de- 
tailed experimental study of Cavaille and Dreschner [Q; 
(d) It is a standard experimental practice to report the 
magnitude of the ion-enhanced diffusion using an effec- 
tive temperature T e ^ at which the substrate needs to 
be heated to obtain the same mobility as with ion bom- 
bardment 0J^] . We can calculate T e ^ using the relation 
D 1 + D exp{-E a /k B T) = D exp{-E a /k B T e ^), that 
has two important consequences. First, the anisotropic 
diffusion constant translates into an anisotropic T e ^ , i.e. 
we have T^f ^ T^' . The experimental methods used to 
estimate T e ^ could not distinguish T^ff and T^f 
However, current observational techniques should be able 
to detect the difference between the two directions. Sec- 
ond, while it is generally believed that ion bombardment 
can only raise the effective temperature since it trans- 
fers energy to the surface, the negative D 1 indicates that 
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along the x direction one could have T e // < T. (e) Fi- 
nally, the results (§)-(U) are based on Sigmund's theory 
of sputtering || that describes sputtering in the linear 
cascade regime. The energy range when this approach is 
applicable lies between 0.5 KeV and 1 MeV, the precise 
lower and upper limits being material dependent. Thus, 
we do not expect (^)-@ to apply to low energy (few eV) 
ion-enhanced epitaxy. 

Ripple formation — Sputtering may lead to the devel- 
opment of a ripple morphology on the surface (l3j |. The 
origin of the ripple formation is an ion-induced instabil- 
ity: valleys are eroded faster than crests, expressed by 
negative v x and v y coefficients in (||) |lCj . At short wave- 
length this instability is balanced by surface diffusion. A 
linear stability analysis predicts that the observable rip- 
ple wavelength is determined by the most unstable mode, 
and has the wavelength £ — 2tt yj D /\v\, where v is the 
largest in absolute value of the negative surface tension 
coefficients. Accordingly, the wave vector of the ripples 
is parallel to the x axis for small 9 and perpendicular to 
it for large 6 Q . The large length scale behavior is de- 
scribed by the noisy anisotropic Kuramoto-Sivashinsky 
equation pjjjl^l , leading to roughenin g |]l5f or the devel- 
opment of coarsening ripple domains pM- 

Quantitative comparison with experiments — At finite 
temperature the total diffusion constant given by D = 
D 1 + Dt- As T decreases there is a critical temperature, 
T c , at which D 1 = D T , so that for T < T c the diffu- 
sion is dominated by its ion-induced component, which 
is independent of temperature, in agreement with the ex- 
perimental results of MacLaren et al. || . Unfortunately, 
for most materials the quantities entering in F, a a and 
Dq are either unknown, or only their order of magnitude 
can be estimated. However, using the expression for £ we 
can express T c in terms of measurable quantities free of 
these constants 



l-{2k B T /E a )ln(£ Ion /e To ) 



where £t is the experimentally measured ripple wave- 
length at any temperature To > T c ; £j on is the ripple 
wavelength in the low temperature regime, T < T c , where 
ion induced diffusion dominates, and therefore £j on is 
independent of T; E a is provided by the slope of ln(^) 
versus 1/T in the high temperature regime (T >> T c ). 
Consequently, all quantities in (||) can be obtained from 
a plot of the ripple wavelength as a function of temper- 
ature, so that (|6|) gives T c without any free parameters. 
Such a plot is provided by MacLaren et al. ||, leading 
to E a = 0.51eV, li on = 0.8^m. Using £ To = 2[xm for 
T = 368K |], we obtain T c = 57°C, which is in good 
agreement with the experiments, that provide T c between 
45 and 60°C @. 
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FIG. 3. Phase diagram for the isotropic case a = fi = 1. 
Region I: v x < 0, v y < 0, D' x > 0, D y > and l x > £ y ; 
Region II: v x < 0, v y < 0, D* x > 0, D y > 0, and l x < £ y ; 
Region III: u x < 0, v y < 0, D' x < and D y > 0; Region IV: 
Vx > 0, v y < 0, Dl < and D y > 0. Note that the phase 
diagram is independent of the precise values of J,p, while the 
e dependence is contained in a a . 

Morphological phase diagram — The detailed morpho- 
logical phase diagram is rather complex if the diffu- 
sion is not thermally activated but ion-induced. At 
low temperatures, when Dt is negligible compared to 
D 1 , the ripple wavelengths are £ ! x = 2-K\J 'D^/\i/ x \ and 

£\ = 2TT^D^/\iy y \. In the following we discuss the de- 
pendence of the surface morphologies on the experimen- 
tal parameters and a ai based on the phase diagram 
shown in Fig. |^. 

Region I — The surface tensions, v x and v y , are nega- 
tive while D x and D y are positive, consequently we have 
a superimposed ripple structure along the x and the y di- 
rections. The experimentally observed ripple wavelength 
is the smallest of the two, and since £ x > £ y , the ripple 
wave vector is oriented along the x direction. The lower 
boundary of this region separating it from Region II is 
given by the solution of the £ x = £ y equation. 

Region II — Here the ripple wave vector is oriented 
along the y direction, since £ x < £ y . This region is 
bounded below by the D x — line. At large length scales 
in I and II one expects kinetic roughening described by 
the KPZ equation [|l]Jl|g. 

Region III — In this region D x is negative, while the 
sign of all other coefficients are the same as in I and II. 
Since both surface tension and surface diffusion are desta- 
bilizing along x, every mode is unstable and one expects 
that the KPZ nonlinearity cannot turn on the KS stabi- 
lization the system being unstable at large length 
scales as well. This instability is expected to lead to ex- 
ponential roughening. The lower boundary of this region 
is given by the v x = line. 

Region IV — Here we have v x > 0, v y < Q, D xx < 
and Dy y > 0, i.e. one expects the surface to be period- 
ically modulated in the y direction, leading to a ripple 
structure oriented along the x direction. In the x di- 
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rection we have a reversal of the instability: the short 
length scale instability generated by the negative D x is 
stabilized by the positive surface tension v x . Thus, there 
is no ripple structure along the y direction. Regarding 
the large length scale behavior, along the x direction the 
surface diffusion term is irrelevant compared to the sur- 
face tension, thus one expects KPZ scaling. However, 
along the y direction the KS mechanism is expected to 
act, renormalizing the negative v y to positive values for 
length scales larger than i* leading to a large wavelength 
KPZ behavior. 

If thermal and ion- induced diffusion coexist, the rip- 
ple wavelengths are given by l x — 2tt[(Dt + D x ) ^v^ 1 / 2 
and i y = 2tt[{D t + Z^)/K|] 1/2 . The phase diagram for 
intermediate temperatures can be calculated using the 
total D. In particular for high T, when Dt >> D x and 
Dt » Dy, the phase diagram converges to the one ob- 
tained in Ref. the ripple orientation being controlled 
by the v x = v y line (dotted in Fig. ||). Thus, with in- 
creasing Dt the phase boundary between the regions I 
and II converges to the v x = v y line and the D x = 
boundary separating the regions II and III shifts down- 
wards, eventually disappearing. However, in the interme- 
diate regions new phases with coarsening ripple domains 
appear as the D x and the v x — lines cross each 
other. 

While we limited our discussion to the effect of the 
ion-induced diffusion on the surface morphology, the re- 
sults (Q)-(||) can be used to investigate other phenom- 
ena as well, such as island nucleation. The experimental 
verification of the above results would constitute an im- 
portant step to elucidate the mechanism responsible for 
ion-induced diffusion, with potential applications to ion- 
enhanced epitaxy. 
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